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In this paper, we study the computability of the initial value problem of Combined Korteweg-de Vries 
(CKdV, for short) equation: + uu^ + u^Ux + Mxxi = 0, u{x,0) — ^(x). It is shown that, for any integer 
5^3, the nonlinear solution operator Kr : H^(R) — > C(M,i/*(R)) which maps an initial condition 
data (j) to the solution of the Combined KdV equation can be computed by a Turing machine. 
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1 Introduction 

Differential equations are very popular mathematical models of real world problems. Not every 
differential equation has a well behaved solution. For those equations whose well behaved solutions 
exist, we are interested in how they can be computed. Thus, the computability of the solution operators 
for different types of nonlinear differential equations becomes one of the most exciting topics in effec- 
tive analysis. This answers questions of the type: is it possible to calculate the solutions of some real 
word problems algorithmically? The answers to these questions are unfortunately not always positive. 
However, there are a lot of very interesting equations whose solutions do exist and can be calculated. 
These equations can be called computably solvable equations, in other words, their solution operators 
are computable. This means that, there are Turing machines which can transfer the initial data to the 
solutions of the equation in some particular spaces. For example, Klaus Weihrauch and Ning Zhong Q 
have shown that the initial value problem of Korteweg-de Vries (KdV) equation posed on the real line 
M: Ut + uux + Uxxx = 0, f,;c G M, m(.x,0) = ^{x) has a computable solution operator. 

In this paper, we investigate a variation of Korteweg-de Vries equation: Uf + uUx + u^Ux + Uxxx = 0. 
This is often called Combined Korteweg-de Vries (CKdV) equation. The Combined KdV equation is also 
an important equation which is frequently used as a mathematical model in physics, hydrodynamics, 
biological and chemical fields. We will show that the solution operator of the CKdV equation is also 
computable. This extends the results of 12]^. The proof of the main theorem is given in Section 2. 

2 Main result 

In this section, we use a similar approach as that in Q and retrieve two estimates to prove the main 
result. We use Type 2 theory of effectivity (TTE) as computation model. More relevant details can be 
found in |7]. 
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Korteweg-de Vries Equation 



For rigorous notation we occasionally write u{x,t) := u{t){x), where u{t) : M — C(M;//*(M)). More 
precisely, we are interested in the following initial value problem (IVP, for short) of CKdV equation on 
the real line M, 



Then we consider the solution operator which maps the initial data (jp G H\R) to the solution 
ueC{R;H'{R)),fors>3. 

The following is the equivalent integral equation of the initial value problem ([T]) 



(2) 



where m(0(x) := u{x,t) ,E (t) (x) := e'""' , and ^ {(p) {x) = J e-''^ (p{^)d^. 

R 

We use the following iterative sequence with the initial data <p as the seed: 

r voit) = ,^-'iEit).^{cp)) 

\ V;+i (0 = Vo (0 - /o' - T) • ^ {^vj (T) + iv^ (T)) ) ) dT. 

The iterative sequence ^ is contracting near f = 0, thus the sequence converges to a unique limit. 
Since the limit satisfies the integral equation it is the solution of the initial value problem ([T]) near 
f = 0. To prove that the solution operator is computable, we need to construct a type-2 Turing machine 
to compute it. 

Firstly, we define the operator: 

S{u, (p,t) = ^-\E{t) . ^(<p)) - {E{t-r).^ {£ (^) + 1 „3 (^))))dT 

which is ([p — )• 5s], 5i,p, 4)-computable. This follows from Lemma 3.2 in Q straightforwardly. There- 
fore, the function S{u,(p){t) := S{u,q),t) is ([p —?- 6<,],ds, [p — ^ 5^]) -computable. Then we define the 
function v : S{R) x N ^ C{R : S(R)) by 

v((p,0)=5(0,(p) 
v{(p,j+i)=S{v{(p,j),(p). 

It is easy to verify that v is {5s, Ym, [p — ^ 4])-computable. 

Now we can show several propositions which lead to the proof of our main theorem. 

Proposition 2.1 Ifu{x,t) is the solution of the IVP (E]), then there is a computable function e : N x M x 
M — 7- M which is non-decreasing in the second and third argument such that 

sup II u{x,t) 11,^4(11 ^> h), 
where Cj (r) : = e {s, T, r), s is an integer and s ^3. 

Proof. See [3] □ 
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Proposition 2.2 Let := 5(0, <p), and v^+i := S{vj, <p). // 



a^rV2[8(3 + r)3/2 II CP ||2 +4(3 + r) II <p II,] ^ 1 , 



yJ+^{t)-vj{t) ||,^2-^(3 + r) 



1/2 



9 



Proof. 



According to Proposition 2.1 and lemmas4.8 in [7], if T > ,m € X^, we can obtain that 



/ ||mm;c + m^M;c|| dt ^ a^ri/^ll 
Jo " "•^ 



where = ^-2' ■ T^/^ + 1, = {m G C([0, T] (R) ; (m) < 00)} is aBanach space with the norm 
II m||^j (detail see Definition 4.7 in [7]). 

Then let W(?)<p = ^ / e'^^e'^''9(^)d^. 
■\/2k Jr 

Since v° := 5(0, <p) ,v^+i := 5(v^<p), by Lemma4.9 and4.10in [7], for j^l, 
\\yJ\\^. = \\W{t)(p-^-J^W{t~T)[ivJ-')\dz-^J^Wit 



<(3+r)i/^|| <p||, + (3 + 



V V 



■ dz 



II 3 



^ (3 + r)^/2 II ^ 11^ + (3 + J) 1/2(^.^1/2 II ^j-i ||2^ +(3 + r)i/2(^s ji/2 II 11^^ _ 

Let r > , such that of ri/2[8(3 + r)^/^ || <p ||2 +4(3 + 7) || <p ||,] ^ i. From || ||x^=|| W (?) <p ||x^< 
(3 + r)^/^ II <p II, we obtain by induction 

II v^' ||x^^2(3 + r)i/2 \\(p\l.iforalljeN) 



For ^ 2, 

iiw-v^--i ||^.=||^V(?-T)(i([(v^--i: 



^(3 + r)i/2afri/2 II (v/-i)' + v^'-iv^'-2+(v^--2)2 



-21 



+ (3 + r)'/2^fri/2 II [v^'-i+v''-2] . [vJ-^-vJ 

^ {3 + TfH(4T'^^\\(pf, II yj-'-yj-^ ||z| + (3 + r)4c4ri/2 ||(p||^ || v^'-^-v^'-^ \\x^ 
^ _ II v7 11^,. 

If af ri/2[8(3 + r)3/2 II (p ||2 +4(3 + r) II 9 II,] ^ i, then we obtain the result that 

II v^-+i(0-W(0 11.^11 v^'+i(0-v^'(0 lb<2-^--i(3 + r)i/2 II <p II, . 



□ 
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Korteweg-de Vries Equation 



Proposition 2.3 Letv{t) = S{v,q)){t) , v„ (f) = 5'(v„,(p„) (f) • ^ 

afri/2(3 + r)3/2(i2 II (p ||2 +16 II (p II, +6) ^ ^, 

then we have 

\\v{t)-vn {t)l^2{3 + Ty/^\\(p-(p„l, 

where = ^ -T -T^l^ + \ . 

Proof. Since v{t) = S(y,(p) (t) , v„ (t) = 5(v„,(p„) (?) , by Lemma 4.8 in UJ, we obtain the result as 
following: 

\\v-Vn\\x^ = \\W{t){(p-(p„)-y^W{t-r) [v^-v^^dT-y^W{t-T) [v'-vl]^dr\\x^ 
^ (3 + r)i/2 \\(p-(p„ II, +(3 + r)'/'afr'/2 y ^2 ^^^^^^2 . y v-v„ lUf 
+ (3 + r)^/2afrV2 II v + v„ llx^ • ||v-v„ lUf . 

By Proposition 2.2, if af.r^/2(3 + 7)^/2(12 || (p ||2 +16 || (p ||, +6) ^ 5,notice that || % ||,^|| (p ||, +1, 
then 

II V - V„ llx^ ^ (3 + r)'/2 \\(p-(Pn lis- +12 (3 + r)^/2 7^1/2 II <p ||2 . II V - v„ ||x^ 

+ 4(3 + r)afrV2||,p||^^.||v-v„||xf 

^ (3 + r)^/^ \\(p-(Pr, II, +i II V-V„ llz^ . 
1 /2 

Therefore || v — v„ \\x^^2{3 + T) ' || (p — % ||,, the sequence {v„} is uniform convergence. □ 

Finally, we give the main results are as follows: 

Theorem 2.4 The solution operator : //^ (M) — )■ C(M;//'^(M)) of the initial value problem problem 
^ is {5f{s, [p —7- dH'])-computable for any integer s ^3. 

Proof. For a given initial value cp G H\M.) and a rational number f > we will show how to compute 
the solution u{t) of the initial value problem ^ at the time interval ^T. For this purpose, we first 
find some appropriate rational number T such that < T < f , and show how to compute u{t) from t' 
and Xj/ := u{t') at the time interval [t',t' + T], ^ f' ^ f, by a fixed point iteration. Using this method, 
we can compute the values u {T /2m) successively for m = 1,2, •• • and finally u (t) for any ^ f ^ f. 
If Ut + uux + u^Ux + Uxxx = 0,u{x,t') = \if{x), and v is defined by v{x,t) := u {x,t + t') , then 

( Vt + vVx + v^Vx + Vxxx = xeR,t^O, 
\v{x,0) = Y{x). 

We assume that the initial value y £ H'\W) is given by a 5//!-name, i.e., by a sequence i/Zq, V'^i , • • • 
of Schwartz functions such that || xjf — Vn Il.v^ 2^". For any « G N, we define function v°,v,^,--- in 
C(M : 5(M)) by 

vl:=S{0,Yn), vi+':=S{vlwn). 

We note that the sequence {v/,} can be computed from i//,,. By Proposition 3.3, the iterative sequence 
v°, v,^, • • • converges to some v„, then v„ is the fixed point of the iteration S and satisfies the following 
internal equation: 

V„ {t) =S{Vn,¥n) = ^-\E{t) . ^{(p)) - f ^-^{E{t-T)-^ ii (t) + 1 „3 (^))))dT 

t/ 
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hence v„ solves the mitial value problem: 

dVn , dv„ 2(^^n dv„ 

~dr "~dx '''~dx 'dx^^ v„(x,0) = V/-„(;c). 

By Proposition 2.3, we will show that, by a contraction argument, for some sufficiently small computable 
real number T > (depending only on (p and T), vi (t) — )■ v„ (?) as 7 — oo for all n, and v„ (t) — v (f) as 
n ^ 00^ sufficiently fast and uniformly in ? € [0,T]. We recall that v is the solution of the initial value 
problem (O. Then we can effectively determine a computable subsequence of the double sequencejvii } 
which will converge fast to v uniformly in f G [0,T]. 

Since v is the limit of a fast convergent computable sequence, v itself is computable. So Kr : 
{(p,t) ^ u{t) is (5//i,p,5//s) -computable for f ^ 0, the reflection R : SiM) — )• S{M),R{'^f){x) := ^{—x), 
is (5//i, 5//i)-computable. Define u' [t) (x) := u {—t) {—x). Then u[ + m'm!^ + u'^u'^ + m^-^ = and for 
t ^ 0,m(-0 =/?om' (?) =/?o7S:m(m'(0),0 =RoK^{R{(p),t), i.e., m(0 = Ro K^{R{(p) ,-t) for? ^ 0. 
Therefore, as the two computable functions join at 0, Kg^ is computable for ? G M. (see (71) □ 

Thus, we prove the main result, and we can see that the machine searches for fast approximations 
to u{x;t), and computes the solutions of the Combined KdV equation with arbitrary precision. This 
approach can be extended to other nonlinear equations. 

References 

[1] J. Barros-Neto. An introduction to the theory of distributions. Pure and Applied Mathematics, 14(1993). 

[2] D. Lu and Q. Wang. Computing the Solution of the m-Korweg-de Vries Equation on Turing machines. Theoret. 
Comput. 5c/.,202:219-236(2008). 

[3] D. Lu and R. Zheng. The boundness of the solution of the Combined KdV equation under the Turing com- 
putable circumstance. Mflf/!.Ap/5/.,21(3):814-818(2008). 

[4] K. Weihrauch. Computable Analysis. Springer Verlag, Beriin, Heidelberg 2000. 

[5] N. Zhong and K. Weihrauch. Computatbility theory of generalized functions. J. Assoc. for Computing Ma- 
chinery, 50(4): 469-505 (2003) 

[6] J.L. Bourgain. Fourier transform restriction phenomena for certain lattic subsets and application to non-linear 
evolution equation ,Part II: the KdV equation. Geom. Funct. Ana/., 3:107-156(1993). 

[7] K. Weihrauch and N. Zhong. Computing the solution of the Korteweg-de Vries equation with arbitrary preci- 
sion on Turing machines. Theoret. Comput. Sci., 322:337-366(2005). 

[8] K. Weihrauch and N. Zhong. Is wave propagation computable or can wave computers beat the Turing machine? 
Proc. London Math. 5oc.,85(2):3 12-332(2002). 

[9] K. Weihrauch and N. Zhong. Is the linear Schrodinger Propagation Turing Computable? CCA 2000,369-311 . 



